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1. Introduction. 

Nowadays, it is well established that the chaotic behavior is the intrinsic property of some of the 
nonlinear dynamical systems [1-6]. Due to the nonlinearity most of the rate equations (charge- 
carrier dynamics) in semiconductor physics, such an unpredictable behavior also can be easily 
detected in most of the realistic models of photoconductors. The technological importance of 
photoconductivity, particularly its applications in high-speed and far-infrared photodetectors and 
similar devices, leads to the need for a careful investigation of the charge-carrier dynamics. One of 
such models was proposed recently in [7], see also [8-9]. To be more specific, that model consists 
O^pf three coupled nonlinear ordinary differential equations describing the population densities of 

Q-electrons n, trapped electrons m, and holes p in the following form of rate equations 
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O^vhere: G is the rate of photoexcited charge carriers;ci and C2 are the capture rates of active 
Ofcenters for electrons and holes; N t is the total number of traps in the system; m is the number 
Q^)f occupied traps at time *; ot\ is the probability of capture of the electrons from the conduction 
~^band to the traps; 71 is the probability of ejection of electrons from the traps to the conduction 
^tand by thermal excitation; 8q is the recombination constant for captured electrons. Let 

. n(t — 0) = ni, m{t — 0) = mi,p(t — 0) = pi, (2) 

c : 

Cbe the initial conditions. 

Hn dealing with such systems (1) it is always convenient to use dimensionless systems: dimension- 
• Jess variables, time and constants are introduced by the following relationships: 

X: 

x = nn \ y = mn x , z = pn , t x = tt \ 

_ _ 1 

h = Ci 1 ,g = Gt riQ l , 

n = (aiio) -1 , 
k = ctiNtt ,at = 7i* , 

W = c 2*0) 

[3 = 5 t n 
= 5 ai\ (3) 
The corresponding dimensionless system is of the form: 

_ = g - (k + l)x + xy + ay, 
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dy , , . 

— = kx-xy- fjyz - ay, (4) 

^7 = 9 ~ Pyz - kz, 

According to [7-9], nonlinear system (4) exhibits chaotic behavior with some values of systems 
parameters g,k, [3 , a ,k. Namely, using the data presented in [7]: G = 10 24 cm _3 sec _1 (for 
the correct dimension purposes we use this value of G instead of G — 10 16 cm~ 2 sec _1 as in [7]); 
N t = 510 14 cm- 3 ; 7l = 0.83sec _1 ; a = LSIO' 3 sec' 1 ; c 2 = 1.510 _5 sec _1 ; <*i = 2.510~ 15 cm- 3 sec- 1 ; 
5q = 10~ 15 cm~ 3 sec~ 1 we obtain the following values for t ,rio: t = y^,tiq = 610 n cm -3 and 
for dimensionless control parameter g = 1.110 15 and other parameters for chaotic behavior to 
occur:/? = 0.4, k = 10~ 2 , k = 8.3310 2 , a = 5.510 2 . As the chaotic, unpredictable behavior was 
uncontrolable, ungovernable, and so undesirable in some situations, initially researchers tried to 
avoid chaos and to deal with the systems in the range of the parameters, where chaos was not 
generated. The situation has changed dramatically since the discovery of the possibility of chaos 
control by Pecora and Carrol [10] , Ott, Grebogi and Yorke [11] in 1990. These seminal papers 
[10-11] have induced avalanche of research works in the field of chaos control. Chaos synchroniza- 
tion in dynamical systems is one of methods of controling chaos, see, e.g. [10-17] and references 
therein. The interest to chaos synchronization in part is due to the application of this phenomenen 
in secure communications, in modeling of brain activity and recognition processes, etc [10-17]. Also 
it should be mentioned that this method of chaos control may result in improved performance of 
chaotic systems [10-17]. In our case of charge-carrier dynamical system chaos control could be used 
for the optimization of the performance of photoconductor devices. 

According to [10] synchronization of two systems occurs when the trajectories of one of the systems 
will converge to the same values as the other and they will remain in step with each other. For the 
chaotic systems synchronization is performed by the linking of chaotic systems with a common 
signal or signals (the so-called drivers): suppose that we have a chaotic dynamical system of three 
or more state variables. In the above mentioned way of chaos control one or some of these state 
variables can be used as an input to drive a subsystem consisting of remaining state variables and 
which is a replica of part of the original system.In [10] it has been shown that if the real parts of the 
Lyapunov exponents for the subsystem (below: sub-Lyapunov exponents) are negative then the 
subsystem synchronizes to the chaotic evolution of original system. If the largest sub-Lyapunov 
exponent is not negative, then one can use the nonreplica approach to chaos synchronization [15]. 
To be more specific, one can try to make negative the real parts of the conditional Lyapunov 
exponents of the nonreplica response system. As it has been shown in [15] from the application 
viewpoint nonreplica approach has some advantages over the replica one. The above-mentioned 
chaos synchronization method [10] (replica approach) is applied to different chaotic dynamical 
systems, see, [7-11] and references therein. As it is already underlined recently a new approach- 
nonreplica approach to chaos synchronization is proposed in [15]. A detailed analysis of that pa- 
per shows that for high dimensional systems the calculation of conditional Lyapunov exponents 
in general requires formidable, tedious numerical and analytical calculations. 
This paper is dedicated to the chaos synchronization in the charge- carrier dynamical systems in 
photoconductors within both the replica and nonreplica approaches. It has been shown that by 
using the boundedness of the solutions to the dynamical systems, nonreplica approach to chaos 
synchronization and Routh-Hurwitz criteria it is possible to make negative all the conditional Lya- 
punov exponents without complex numerical and analytical calculations. This is the main feature 
of the paper. 

Thus consider the possibility of chaos synchronization in the dynamical system (4). The system 
has the following steady state (fixed point) solutions: 



x gt = g(l+^K- 1 y' t )- 1 , 



a(3K-\y st ) 2 + (a + g 

+P^ 1 g)y st -kg = 0, (5) 

For the begining we investigate system (1) qualitatively, as it has been shown in [7], for different 
values of parameters the system exhibits different types of solutions, such as oscillatory, chaotic, 
etc. After linearization about the steady state solutions (5) of the system (4) we obtain that 
characteristic equation's roots are satisfying the following equation: 

A 3 + aX 2 + bX + c = 0, (6) 

where 

a = k + f + x st + (0 - l)y st + a + k + f3z s \ 
c={k + (3y st ){a + x st ) 
+K(k - y st + l)(3z st , 
b=(k-y st + l)(K + f3y st + (3z st ) 

+f3KZ St 

+ (a + x st )(l + K + Py st ), (7) 
In order to have oscillating solutions, the value of 

„ c 2 b 3 ca 3 

F = 1 1 

4 27 27 



abc a 2 b 2 



(14) 
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should be positive [18]. The period of the oscillations could be estimated by the following relation- 
ship: 



T = 3- a5 7r((-| + F - 5 ) 



where 



-(- q - 
y 2 

-F°- 5 )i)-\ (15) 

q a 3 ab 
2 ~ 27 ~ Y 
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As the estimations show the period of oscillations approximately conform to the result of [7]. Now 
suppose that the nonlinear system's (4) parameters values are so that system (4) exhibits chaotic 
behavior and we will try to explore the possibility of chaos synchronization in the sense of Pecora 
and Carroll [10]. First consider the x variable as a driver. Then variables y, z form the replica 
response system (with the subscript "r"): 

— — = kx — xy r — (3y r z r 
at i 



-ay r , (9) 
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A 2 + \(x + (3(y + z) + a 



+k) + (K + yp)(x + a) 
+ K (3z = 0, 



(10) 



Here x(t),y(t), z(t) are the solutions of the system (4). 

As it follows from the physical meaning of x, y, z they should be positive and bounded. The 
boundedness of the solutions of the nonlinear system (4) also follows from the Lorenz criteria (for 
more details, see [19]). It means that the roots of equation (10) are negative, and therefore chaos 
synchronization is possible in the case of x driving. Similar calculations show that in the case of 
y, z drivings chaos synchronization also could be realizable. Indeed, for the case of y driving the 
eigenvalues of the Jacobian matrix of the corresponding replica subsystem are calculated from the 
relationship: 



As the calculations show for the above given values of control parameter g and other parameters 
k > y. It means the both A's are negative and therefore chaos synchronization is possible. Taking 
into account the above mentioned inequality k > y one can easily establish that in the case of z 
driving chaos synchronization also is possible, as under this inequality the roots of equation 



are also negative. Thus we have established that the replica approach is sufficient to perform chaos 
synchronization in the case of nonlinear system (4). As it was underlined above, the nonreplica 
approach to chaos synchronization could be used not only in the case of failure of replica approach, 
but also in the case of success by the replica approach to perform chaos synchronization. In the 
latter case one can use the nonreplica approach, say to make larger (in magnitude) the negative 
conditional Lyapunov exponents, as the very important quantity of chaos synchronization time 
inversely proportional to the largest conditional Lyapunov exponent. The application of the 
nonreplica approach (in our case) also could be justified in situations when conditions k > y, k+1 > 
y do not take place. (According to the definition of k (see (3)) such a situation could be realizable 
by changing the values of ai,N t ,ci). As it can be seen easily, then in the case of y driving within 
the replica approach, one of A's become positive, and therefore chaos synchronization is impossible. 
Then to perform chaos synchronization we could recourse to the help of the nonreplica approach. 
As the calculations show under the above mentioned conditions k < y, k + 1 < y the realizability 
of chaos synchronization in the case of z driving is also under question. Here we will restrict 
ourselves to the nonreplica approach in the case of z driving. The simplest nonreplica response 
system [15] could be written in the following manner: 



(\- K -yP)(\ + k + l-y)=0, 



(11) 



A 2 + X(x + /3z + a + k + l-y) + (3z(k + 1 - y) 



+x + a = 



(12) 



dx, 



nr 



g-(k+ l)x nr + x nr y + ay 



dt 



i 



+si(y, 



'nr 



y)=A 



i 



dy nr 
~dU~ 
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-ay + s 2 {y„ r - y) = A 2 , (13) 
■ g f3yz nr Kz nr 
+s 3 {y nr -y)= A 3 , 

where s±, s 2 , s 3 are some arbitrary parameters (constants). The eigenvalues of the Jacobian 

<"» 

(where A = (A±, A 2 , A 3 ), x = (x nr , y nr , z nr )) are to be found from the equation: 

A 3 + \ 2 {k + l-y + K + Py-s 2 ) 

+X(k + l- y )( K + Py-s 2 )- s 2 (k + fly) 

si(k - y) + f3ys 3 ) 

+(k + 1 - y)((3ys 3 - s 2 (k + (3y) - si(k 

-y)(K + Py) = 0, (15) 

Our aim is to make negative the roots of this equations without explicit calculations. To achieve 
our goal we can use the following advantages: 1) As we have mentioned above the solutions of the 
nonlinear system is bounded; 2) the equation for A contains some arbitrary parameters s±, s 2 , s 3 . 
In addition, there is no need to calculate these roots explicitly (in the case of higher dimensional 
systems the task is rather formidable) and investigate their dependence on the arbitrary constants 
with the aim of making the roots negative. For this purpose one can just use the Routh-Hurwitz 
criteria: the necessary and sufficient conditions for the equation (15) to have roots with negative 
real parts are : 

d\ > 0, a 3 > 0, a\a 2 — a 3 > 0, (16) 

where ai,a 2 ,a 3 are coefficients of A 2 , A 1 , A in the equation (15), respectively. If the nonlinear 
chaotic system's parameters value permits the condition k > y, then as it is clear from the equa- 
tion for A , just the presence of larger negative s 2 is sufficient to make negative the roots of the 
equation (15). In other words, we can safely put Si = s 3 = O.If k — y < Q,k + 1 — y < Q then 
larger negative s 2 and larger (in comparison with the absolute value of s 2 ) positive s\ values with 
S3 = could make realizable chaos synchronization within the nonreplica approach. Without 
doubt, the presence of all the arbitrary constants si,2,3 in principle further alleviates the task of 
making negative the roots of the equation (15). 

Thus in this paper we have shown that using the boundedness of the solutions of the some of 
the charge-carrier dynamical systems, nonreplica approach to chaos synchronization and Routh- 
Hurwitz criteria one can make negative the real parts of the conditional Lyapunov exponents 
without formidable, cubmersome and tedious numerical and analytical calculations. (For more de- 
tailed presentation of this approach to chaos control, see [20]). 
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